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NUMERICAL SOLUTION OF HYPERBOLIC EQUATIONS 
AND SYSTEMS BY A METHOD OF THE 
RUNGE-KUTTA TYPE. II 

by 

Nguyen Kong Tuy 


In the first part of the article the two-iteration algorithms of the Runge- 
Kutta type were applied to the solution of the Cauchy problem for hyperbolic 
equations and systems with two independent variables , where the initial data 
are given along the line segment x + y = const. 

In the second part of this article an analogous problem is discussed for 
one equation with Cauchy data along the curve segment. 


Below are used designations of the article [1] and the numbering of its 
paragraphs, formulas and figures is continued. 


6. Statement of the Problem 

Let there be an equation 

= f(x,y,u,p,q) (41) 

^u is an unknown function of x,y; p = u^, q = u^with initial conditions given 
along the segment AB of a certain curve T 0 in the form of: 

u° = u°(x), p° = p°(x), q° = q°(x) ; (42) 

u °> P°, q°, f are assumed to be continuous and differentiable a sufficient number 
of times. The term T 0 is assumed to be a sufficiently smooth line given in the 
form 


y = y°(x) or x = x°(y) 
while x°(y) , y°(x) are decreasing functions. 


(43) 
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Let us construct in the region of ABC of the definition of the solution a 
grid in the following (Figure 3) . Let us subdivide AC into sufficiently small 
intervals of length h; from the points of division we will plot straight lines 
parallel to the axis Oy, and from their intersections with r 0 , those parallel to 
the axis Ox, and we will distribute the nodes of the grid according to layers. We 

th 

will connect the nodes of the n layer by the curve r . The curve r can be 

n n 

considered a segment of the initial curve r 0 displaced by nh in the direction Ox 
so that we have the equation T in the form 


y n (x) = y°(x - nh) = y°(z) , z = x - nh . 

til 

In making calculations of the n layer we will use, for simplicity, the equation 

of the curve r - in the form (43) , and we will denote the earlier found values 
n - 1 

st 

u, p, q in the nodes of the (n - 1) layer, as well as the initial data, by 
„o o o 
u p , q . 

For an elementary curvilinear triangle MNP of the grid with vertices 

M(x 1 ,y 1 ), N(x 2 = Xj + h, y 2 = y°(x 1 + h) ) on the (n - l) St layer and the 

vertex P(x 2 ,y 1 ) on the n^ 1 layer (Figure 4) , the following relationships are 
known: 

Xj + h 

u(P) = u°(Xj) + J p°(x)dx + JJ F(x,y)dxdy , (44. a) 

y°(x!) 

p(P) = p°(x 2 ) + / F(x, + h,y) dy , (44. b) 

y°(x t + h) 

Xj + h 

q(P) = q° ( Xj ) + J F(x,y 1 )dx , (44. c) 


where the multiple integral in (44. a) is taken over the region MNP. As in [1] , 
for the increments A* u = (P) - u(M) , A*p = p(P) - p(M) , A* q = q(P) 

- q(M), by using (44) we can find expansions in powers of h at the point M: 
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y2 = y°(* i *»)] 

FIGURE 4 


A*u = hp°(x 1 ) + ~ ^ 

h 2 dy°(x 1 ) 


h 2 dp°(Xj) h 8 d 2 p° ( Xj ) 


31 dx 2 
h 3 d 2 y°(x 1 ) 


2-. dx F(M) -fr :s i5r-F(M) 
h 3 dy° ( Xj ) 2 . 8F(M) , dy°(xi) 9F(M) 


31 dx 


ax 


dx 


3y 


+ o(h 4 ) , 
(45. a) 


, . dp°(x 1 ) h 2 d 2 p°(x 1 ) dy°(x 1 ) 

A *^ !Sh -^r L+ -2i dbT 2 - h -^ r ' F(M) 


h 2 d 2 y° ( Xj ) 
21 dx 2 


F(M) 


dy°( Xl ) 
21 dx 


/ 2 aF(M) + dj^9F(M)\ + 3) f 

y ax dx ay J ; 

A *q = hF(M) +-|7^^ + 0(h 3 ) . 


(45. b) 
(45. c) 


We set forth the problem of the construction of algorithms of the Runge- 
Kutta type which give approximate values of the transformations Au, Ap, Aq, 


3 


7„ Formal Algorithms of the Runge-Kutta Type 


x l " h 

We will denote j p°(x)dx = 1(h). (For calculations of 1(h) 
X) 

see Paragraph 4) . We introduce values of k: 


k = k(h) = y°(x 1 ) - y Q (xj + h) (h > 0) . 


Upon the assumption that y°(x) is a decreasing function, the values of k 
will be positive. In Figure 4, k(h) is equal to the length of NP. (Construction 
of the grid in Paragraph 6 and introduction of values of k(h) is borrowed 
from [2], From the determination of k(h) we have: 


k(h) = 0(h) 


k'(0) 


dk(h) 

dh 


h = 0 


dy°(xi) 

dx 


Assume 


k" 


( 0 ) = 


d 2 y 0 (Xi) 

dx 2 


hk_ + h 3 dVtxQ . 
2 12 dx 2 


(46) 


Let us note that S is equal to the area of the curvilinear triangle MNP, 
correct to the order 0(h 4 ) . 

At the first iteration we assume that: 

u 1 - u°(x,) + 1(h) + SF, , 

Fj = cojF(M) + tjF(N) , co j , r j — 0 , coi + t \ 

p 1 = p°(x 2 ) + kFf , F^ - F(N) , 

q 1 = q°(x,) + hFf , Fj = F ( M) . 


= 1 , (47. a) 

(47. b) 
(47. c) 


We form the elementary increments: 

A*u = u 1 - u°(x,) > 
A J p = p 1 - p 0 (Xj) > 
A J q = q 1 - q°(Xi) ; 


(47. d) 



(47. e) 


5 ] u = u 1 - u°(x 2 ) > 

<5 ! P = P 1 - P°(x 2 ) , . 

5^ - q 1 - q°(x 2 ) 

Furthermore, we introduce the intermediate points M 2 (j = u,q) , 
No(j = u,p): 


M 2 = ^Xj + p J h ; yj + cr J k ; u^Xj) + o^u ; 

p°(x,) + /J-’a'p ; q°(xj) + 7 ^^) , (48. a) 

N 2 = (x 2 + p J h ; y 2 + ct ‘'k ; u°(x 2 ) + cv S*u ; 

p°(x 2 ) + /3 Vp ; q°(x 2 ) + 7~ ;i 6 1 q) , (48. b) 

where p\ a\ a ^ , . . . , 7 •* are numerical parameters. 

At the second iteration we assume that: 

u 2 = u°(Xj) + 1(h) + SF 2 , 


F, = 


*)- 

► t 2 f(n 2 ) 


W 2 ,T 2 ^ 1 

0 , U) 2 + T 2 = 1 , 

(49. a) 

P 2 = 

: P°(x 2 ) + 

kF 2 

F 2 

= f(n?) , 

(49. b) 

q 2 = 

= q°(X!) + 

hF 2 " 

f 2 = 

f(m?) , 

(49. c) 


A 2 u = 

u 2 - 

U 0 ^) 




A 2 p = 

P 2 - 

P # (Xi) , ' 

► • 

(49. d) 
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a* 

CM 

<1 

q 2 - 

q°(x,) 




Finally, for Au, Ap, Aq we take the following linear combinations: 

Au = XjA^ + A 2 A 2 u ,' 

Ap = PiA J p + M2A 2 p , > 

A p - A 1 q + y 2 A 2 q 


( 60 ) 
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(V <V v i 


are numerical parameters ) . 


8. Numerical Determination of Parameters 


For A 1 u, A 1 p, A*q are valid the following expansions in powers of 
h at the point M: 

A‘u - h P «,x,) + | + | iQp ♦ £ k.(0)F(M, 


21 


_h 3 

31 


| k " (0)+ i^^al F(M) 


3 / dF.\ 

+ -^- 3k'(0) (-^) + 0(h 4 ) , 

3. \ an / M 

A*p = h d2 ^ Xl) + hk* (0) F(M) 


(51. a) 


+ Y k"(0)F(M) +^r 2k'(0)^pj + 0(h 3 ) , (51. b) 


dF? 


h 2 / dF “\ 

A*q = hF(M) +— 2\-T~j + 0(h 3 ) . 

3- \ an / M 


(51. c) 


By comparing (50) and (51) and substituting k' ( 0) , k"(0) in (46) we 
obtain expansions in powers of h at the point M for Au, Ap, Aq: 


Au = (\i + X 2 ) 


o,_ x ^ A 2 d P°( x i) + h 3 d 2 p°(xi) 


hp (Xl) + ~2[ —&T- ' 3l dx2 


h 2 dy°( Xl ) 


- (Xl + X2) ii ri dr iiF(M) - (Xl + Xz) T: dx? 

h 3 dy°(x 1 ) _ /. dFi . dF ‘ 


h 3 d 2 y°(x 1 ) 


F(M) 


or ^ x 3 f x ! + x 2 “Sr ) + °(H 4 ) > (52. a) 

31 dx l dh 2 dh M 
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Ap = (Mi + M2) 
- (Mi + M2) 


h dpVi) , h 2 d 2 p°(x 1 ) 
dx 21 dx 2 


h ^r^ F(M) "li F(M > 


dx 2 


,2 1 / d F ^ dF; 


h dy u ( Xl ) 2 


2'. dx 


Ml ir + M2 + 0(h3) • (52,b) 


Aq = (i>i + v 2 )hF (M) + — 2 Ivj 


dF“ 

dh 


^2 


M 

d lL\ 

dh 


+ 0(h 3 ) . (52. c) 


M 


The selection of parameters is effected so that the corresponding expan- 
sions of (52) and (45) coincide for an arbitrary function F(x,y) and arbitrary 
step h. As a result of comparison of (52) and (45) we obtain the following 
conditions: 


Xj + \ 2 = Mi + M2 = *T + ^2 = 1 > 


3 ( X 1 + X 


dh 


dF 2 

dh 


__ 0 9F(M) ^ dy°(x 1 ) 3F(M) 

— A ■ ' 


'M 


ax 


dx 


9y 


„ / dFr dF 2 \ „ 8F (M) , dy°(x 1 ) 3F(M) 

2 ("■ _ dT + M! "dT/ M = “ST W 


(53) 

(54) 


„ / dF" 

2 h-dF + 


dFf 

dh 


M 


aF(M) 

ax 


(56) 


For the right sides of (54) through (56) we have: 


JLE-iLL+li , 9f 9p ! 9f f 
a x - ax + au p ap ax aq 


_aF_ = _af + _af + 9f_ + a_f 3 cl , 

ay ay + au q ap aq ay 


where 


*V= lim 
3x , 


A* p = dp°(x) _ dy°(x) 
h dx dx 


(57) 
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is 

ay 


lim 
h— 0 


q (Xi) + A*q - q°(X! + h) 


= lim 


1 q°(xi) + A*q - q°(x 1 + h) 


dq°(x) 


dx 


- f 


k/h 

/dy°( x) N 


dx 


From (57) it follows that: 


9F(M) + dy°( Xl ) 9F(M) 
“ax dx ay 


r 2 

af 

, i 

dy° 

af 

L 2 

ax 

dx 

ay 

dy° 

n° 1 + 

af 

L 

dx 

q ; 

ap 

V 

af / 

dq + 

,)] 


aq \ 

dx 

/ 

M 


9f 


2p° 


dp 0 

dx 


dy° 


dx 


( 57 ) 


(58) 


For the calculation of the left sides of (54) through (56) it is necessary 
to find (d/dh) M for F(N) , f(m ] 2 ,)f(n J 2 ): 


dF(N) _ dF( Xl + h, yj - k(h)) 9F(N) dk(h) 9F(N) 

dh 


so that on the strength of (46) and (57) 


dF(N) 
dh 

| a f dp° | a f dq° 
3p dx 3q dx 


af + dy° af | a: 
ax + dx ay ai 


Furthermore, from (51) it follows that: 


(dAV 

\ dh , 

/ dA*q ^ 


= P° (Xj) 


M 


dA*p 

^"dh", 


= FfM) 


and from comparisons of 6^, 6^, 5*q wit 


© 

= -q°. 

M 

/dsV) 

= _ dj 

l dh / 

M 

dS*q\ 

dq°(x t ) 


dh 


M 


dx 


+ F(M) 


From (48. a) , using (60) , we have 


dF^) 


dh 


J M L 


j a f j dy° a f J o 9 f 

> — <t —+ a p — 

ax dx ay au 


/v . d/ A af + j 

+ p \ dx dx / 3p 3q 


M 


(61) 


(62) 


dF^) 

dh 


= Xj + 

h, 

43 

M 

1 

£ 

II 

e-j 

>> 

using (61) , 


u 


J\ 

dy° 3f 

M 

[_l 

p iax l 1 

^ ) 

dx 3y 


/„ o 

, o dy° -j o 

dy° 

\af 

+ 

Ip 

+ q* 5 — - a q 

dx 

i au 


+ (v.'-iiitW 

\ dx P dx / 9p 


53! . y -i 5a! + y-if ) 

1 dx Y dx Y Id q 


(63) 


M 
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Let us consider the condition (54) . Its right side is given in (58) and the 
left side can be expressed in the following form, using for this (47. a) , (49. a) , 
(59) , (62) , and (63): 


3 9x ^ lTl + ^ 2T2 + ^ 2Ct, 2 p + P) 


9 f _ dy u 

+ 3 — (XjTj + X 2 t 2 - X 2 u>2<r - X 2 r 2 a) ^ 

3 f r 

+ 3 — I (Xjn + X 2 t 2 + X 2 oj 2 o;)p 0 

+ (X lTl + X 2 t 2 - X 2 t 2 c 5 ) q° 

„ 9 f r dp 0 

+ 3 9p L lTl + XzT2 + X 2 W2/3) "dT 

- dv° 1 

- (X 2 o> 2 /3 + X 2 t^) fl 

+ 3 Iq [ (XlTl + XzT2 “ X2TiY) "S" 

+ (X 2 co 2 y + X 2 T 2 y)fJ . (64) 

Here, all values are given at the point M and, for simplicity, all indices j = u 
are omitted. 

As in [1] we assume: 

a = a , /3 = /3, y = y , p = p - 1 , cr = a + 1 . (65) 

Taking into account (53) , (65) and equating (64) to the right part of (58) , 
we will obtain for the parameters at Au a system of equations (24) from [1]. 
Thus, the results of Table 1 of the article [1] remain valid for the parameters 
at Au. 


By considering analogously the conditions (55) , (56) we will obtain for 
the parameters at Ap, Aq the systems of equations (25) , (26) from [1], so 
that for the parameters at Ap, Aq the results of Table 2 of the article [1] 
remain valid. 
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The scheme of calculations according to the proposed algorithms will 
have a form analogous to (27) , (28) from [1]. 
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